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Abstract. We consider the Glauber dynamics for the 2D Ising model in a box of 
side L, at inverse temperature j3 and random boundary conditions r whose distri- 
bution P either stochastically dominates the extremal plus phase (hence the quo- 
tation marks in the title) or is stochastically dominated by the extremal minus 
phase. A particular case is when P is concentrated on the homogeneous configu- 
ration identically equal to + (equal to — ). For /3 large enough we show that for any 
e > there exists c = c(/3, e) such that the corresponding mixing time T m i x satisfies 
limL^oo P (Tmix > exp(cl/ e )) = 0. In the non-random case r = + (or r = — ), this 
implies that T m i x < exp(cL e ). The same bound holds when the boundary conditions 
are all + on three sides and all — on the remaining one. The result, although still very 
far from the expected Lifshitz behavior T m i x = 0(L 2 ), considerably improves upon 
the previous known estimates of the form T m i x < exp(cL 3+e ). The techniques are 
based on induction over length scales, combined with a judicious use of the so-called 
"censoring inequality" of Y. Peres and P. Winkler, which in a sense allows us to guide 
the dynamics to its equilibrium measure. 
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1. Introduction, model and main results 

Glauber dynamics for classical spin systems has been extensively studied in the last 
fifteen years from various perspectives and across different areas like mathematical 
physics, probability theory and theoretical computer science. A variety of techniques 
have been introduced in order to analyze, on an increasing level of sophistication, the 
typical time scales of the relaxation process to the reversible Gibbs measure (see e.g. 
|17[ ITi] and the recent work on the cutoff phenomenon for the mean field Ising model 
[15]). These techniques have in general proved to be quite successful in the so-called 
one-phase region, corresponding to the case where the system has a unique Gibbs state. 
When instead the thermodynamic parameters of the system correspond to a point 
in the phase coexistence region, a whole class of new dynamical phenomena appear 
(coarsening, phase nucleation, motion of interfaces between different phases,...) whose 
mathematical analysis at a microscopic level is still quite far from being completed. 

A good instance of the latter situation is represented by the Glauber dynamics for 
the usual ±1 Ising model at low temperature in the absence of an external magnetic 



field (see Section 1.2). When the system is analyzed in a finite box of side L of the 
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d-dimensional lattice Z d with free boundary conditions, the relaxation to the Gibbs 
reversible measure occurs on a time scale exponentially large in the surface L d ~ x |27|l26] 
because of the energy barrier between the two stable phases of the system (see Section 



1.3 for a more quantitative statement). When instead one of the two phases is selected 



by homogeneous boundary conditions, e.g. all pluses, then equilibration is believed 
to be much faster and it should occur on a polynomial (in L) time scale because of 
the shrinking of the big droplets of the opposite phase via motion by mean curvature 
under the influence of the boundary conditions. Unfortunately, establishing the above 
polynomial law in Z rf remains a kind of holy grail for the subject and the existing 
bounds of the form exp(cy / L log(L)) in d = 2 [HI [12] and exp(cL d ~ 2 log(L) 2 ) in d > 3 
|25] are very far from it. 

It is worth mentioning that, always for the low-temperature Ising model but with the 
underlying graph G different from Z d , it has been possible to carry out a quite detailed 
mathematical analysis. The first example is represented by the regular <i-ary tree |18] 
and the second one by certain hyperbolic graphs [5]. In both cases one can show for 
example that the relaxation time or inverse spectral gap of the Glauber dynamics in a 
finite ball with all plus boundary conditions is uniformly bounded from above in the 
radius of the ball, a phenomenon that is believed to occur also in Z rf in large enough 
(> 4?) dimension d. 

Moreover polynomial bounds on the mixing time, sometimes with optimal results, 
have been proved for some simplified models of the random evolution of the phase 
separation line between the plus and minus phase for the two-dimensional Ising model 
(see for instance [7] and [IH]). The latter contribution, in particular, partly triggered the 
present work. There, in fact, the opportunities offered by the so-called Peres- Winkler 
censoring inequality |22 j have been detailed in the very concrete and non-trivial case 
of the so-called Solid- on- Solid model. 

Roughly speaking the censoring inequality (see Section 2.4) says that, when consider- 
ing the Glauber dynamics for a monotone system like the Ising model on a finite graph 
and under certain conditions on the initial distribution, switching off (i.e., censoring) 
the spin flips in some part of the graph and for a certain amount of time can only 
increase the variation distance between the distribution of the chain at the final time 
T and the equilibrium Gibbs measure. Therefore, if the censored dynamics is close to 
equilibrium at a certain time T, the same holds for the true {i.e. uncensored) one. 

The fact that the choice of where and when to implement the censoring is completely 
arbitrary (provided that it is independent of the actual evolution of the chain) offers the 
possibility of (sort of) guiding the dynamics towards the stationary distribution through 
a sequence of local equilibrations in suitably chosen subsets of the graph. Of course the 
local equilibrium in each of the sub-graphs is conditioned to the random configuration 
reached by the dynamics outside it and therefore one is naturally led to consider the 
Ising model with random boundary conditions, a quite delicate topic because of the 
extreme sensitivity of the relaxation or mixing time to boundary conditions (see [H [21 
[31 S] for several results in this direction, some of them quite surprising at first sight). 
Moreover it should also be clear that, in order for the guidance process to be successful, 
the distribution of the random boundary conditions at each stage of the censoring 
should be close to that provided by the stationary Gibbs distribution, a requirement 
that puts quite severe restrictions on the choice of the censoring scheduling. 
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The main contribution of this paper is a detailed implementation of this program 
for the two-dimensional, low-temperature, Ising model in a finite box with either ho- 
mogeneous, i.e. all plus (all minus), boundary conditions or, more generally, random 
boundary conditions that are stochastically larger (stochastically smaller) than those 
distributed according to the plus (minus) phase. 

In order to state precisely our results we need to define the model, fix some useful 
notation and recall some basic facts about the Ising model below the critical tempera- 
ture. 

1.1. The standard Ising model. Let A be a generic finite subset of 1? . Each site x 
in A indexes a spin o~ x which takes values ±1. The spin configurations {cr x }xeA have a 
statistical weight determined by the Hamiltonian 

H T (a) = -- Yj °x°y- Yl axT v> 

\i— i/|=i |«— a|=i 

where r = {T y } y £A c are boundary conditions outside A. 

The Gibbs measure associated to the spin system with boundary conditions r is 

y a e n A := {-1, +1} A , 7T T A (a) = exp (-/?#») , 

y3,A 

where (3 is the inverse of the temperature (/? = ^) and A is the partition function. If 
the boundary conditions are uniformly equal to +1 (resp. —1), then the Gibbs measure 
will be denoted by -k\ (resp. n A ). If instead the boundary conditions are free (i.e. 
T y = Vy) then the Gibbs measure will be denoted by 7r£. 

Remark 1.1. Sometimes we will drop the superscript r and the subscript A from the 
notation of the Gibbs measure. 

It is useful to recall a monotonicity property of the Gibbs measure that will play a 
key role in our analysis. One introduces a partial order on Q\ by saying that a < r\ if 
o~x < f]x for all x G A. A function / : i— ► K is called monotone increasing (decreasing) 
if o < i] implies f(a) < f(rf) (f(cr) > f(i]))- An event is called increasing (decreasing) 
if its characteristic function is increasing (decreasing). Given two probability measures 
fi, v on Qa we write \x ^ v if fJ>(f) < for all increasing functions / (with fx(f) we 
denote the expectation of / with respect to fi). In the following we will take advantage 
of the FKG inequalities [TT] which state that 

• if r < t', then tt\ X tt\ 

• if / and 5 are increasing then ir\(fg) > n\{f)n\(cj). 

The phase transition regime occurs at low temperature and it is characterized by 
spontaneous magnetization in the thermodynamic limit. There is a critical value (5 C 
such that 

V/3 > Pc, lira 7rf (cr ) = - lim vr7(cr ) = ma > . (1.1) 

A^Z 2 A A^Z 2 A 

Furthermore, in the thermodynamic limit the measures 71"^" and converge (weakly) 
to two distinct Gibbs measures 7r^ and 71"^ which are measures on the space Z 2 = 
{— 1, +1} Z . Each of these measures represents a pure state. 
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The next step is to quantify the coexistence of the two pure states defined above. 
Let Ax, = { — LL/2J , . . . , [L/2\ } 2 , let ft be a vector in the unit circle S and 0^ the angle 
it forms with e\ = (1,0) and finally let r be the following mixed boundary conditions 



if ft ■ y ^ 0, 
if n • y < 0. 

7± 



The partition function with mixed boundary conditions is denoted by L (n) and the 
one with boundary conditions uniformly equal to +1 by Z^ L . 

Definition 1.2. The surface tension in the direction orthogonal to ft £ § is an even 
and periodic function of of period tt/2, and for — 7r/4 < < 7r/4 it is defined by 

f-\ r cos (0n) , q n 

r,(n) = ton -- — log . (1.2) 



We refer to [21] for a general derivation of the thermodynamic limit (1.2). With this 
definition, one result (among many others) concerning the coexistence of the two phases 
can be formulated as follows [23, . Let m\ L (a) = Y^xeA L a % be * ne total magnetization 
in the box Ax. Then 



ilog(7r{ L (Lm Ai /2j =0)) =r /3 (1.3) 
where rg is the surface tension in the horizontal direction e\. 



lim 

L— >oo L 



1.2. The Glauber dynamics. The stochastic dynamics we want to study, sometimes 
referred to as the heat-bath dynamics, is a continuous time Markov chain on Q\, re- 
versible w.r.t. the measure 7r^, that can be described as follows. With rate one and for 
each vertex x, the spin a x is refreshed by sampling a new value from the set {—1, +1} 
according to the conditional Gibbs measure ir x := 7T^(- | a y , y ^ x). It is easy to check 
that the heat-bath chain is characterized by the generator 

(a/)w = Ew/)-/wi (i- 4 ) 

where TT x (f) denotes the average of / with respect to the conditional Gibbs measure 
ir x , which acts only on the variable a x . The Dirichlet form associated to L\ takes the 
form 

£KfJ) = }Z^U^(f)) 

where Var x (/) denotes the variance with respect to tt x . 

We will always denote by \x% the distribution of the chain at time t when the starting 
point is a. If a is either identically equal to +1 or —1 then we simply write [if or 
fj,^. The boundary conditions r are usually not explicitly spelled out for lightness of 
notation. Sometimes we write ^ t wri en we wish to emphasize that we are looking at 
the evolution for a system enclosed in the domain A. 

The Glauber dynamics with the heat-bath updating rule satisfies a particularly useful 
monotonicity property. It is possible to construct on the same probability space (the 
one built from the independent Poisson clocks attached to each vertex and from the 
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independent coin tosses associated to each ring) a Markov chain {r]t' T }t>0i r) G 
Q\ x riA c , such that 

• for each r G £l\c and cr G f2 A the coordinate process {r^' T )t>o is a version of 
the Glauber chain started from a with boundary conditions r; 

• for any t > 0, % ' < 7y t ,r whenever a < a' and r < t'. 

It is possible to extend the above definition of the generator £ A directly to the whole 
lattice Z 2 and get a well defined Markov process on 0,%? (see e.g. [13 ). The latter will 
be referred to as the infinite volume Glauber dynamics, with generator denoted by C. 

Two key quantities measure the speed of relaxation to equilibrium of the Glauber 
dynamics. The first one is the relaxation time ^reiax* 

Definition 1.3. T re i ax is the best constant C in the Poincare inequality 

Var A (/) := Var ff r (/) < C£ T A (f, f), V / : ft A ^ R. (1.5) 

In particular, for any / : £l\ i— ► R, it follows that 

Var A (e* £ A/) 1/2 < e -*/T rcla x Va r A (/) 1/2 . (1.6) 

We will write gap := gap A for the inverse of T re i ax . 

Another relevant quantity is the mixing time which is defined as follows. Recall that 
the total variation distance between two measures fx, v on a finite probability space Q 
is defined as 

IIm-HI-^EI^)-^)!- M 

Definition 1.4. For ony e G (0,1), we define 

T mix (e) := inf{t > : sup - vr A || < e}. (1.8) 

a 

When e = l/(2e) we will simply write T m j x . 
With this definition it follows in particular that (see e.g. |14j ) 

sup - ttXH < (2e) L ' /Tmix{e)J Vt > 0. (1.9) 

As it is well known (see e.g. (14] ) the following bounds between T re i ax and T m ix hold: 

Frelax < F mix < log (^ T re i ax (1.10) 

where tt* = min CT 7r A (er). Notice that 7r* > e _c ' A ' for some constant c = c(/3) and 
therefore the two quantities differ at most by const x volume. 
Another definition we will often need is the following: 

Definition 1.5. Let [/,, v be measures on J7 A; let a G and V C A. Then, — v\\y 
denotes the variation distance between the marginals of \x and v on £ly, and ay the 
restriction of a to V . 
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1.3. Main results. Our main result considerably improves upon the existing upper 
bound on the mixing time (and therefore also on the relaxation time) when A is a 
square box and the boundary conditions r are homogeneous i.e. either all plus or all 
minus. As a by-product we also get a new bound on the time auto-correlation function 
of, e.g., the spin at the origin for the infinite volume Glauber dynamics started from 
the plus phase 7r+ . Before stating the results we quickly review what was known so 
far. In what follows A^ will always be a L x L box. 

When the boundary conditions are free, a simple bottleneck argument proves that 

T 



T re iax> £2 [f AL (lm AL /2\=0) 



so that (recall (1.3)) 

lim - log(r re i ax ) > Tp. 

L— >oo Li 

In [16] such a result was improved to an equality for large enough values of (3 and in 
for any (3 >(3 C . 

Quite different is the situation for homogeneous boundary conditions, e.g. all plus, 
for which the bottleneck between the two phases is removed by the boundary conditions 
and the relaxation process should occur on a much shorter time scale. In this case one 
expects a polynomial growth of both T ro i ax and T m [ x of the form 

^relax ~ T m i x ~ L . 

The reason behind the difference in the power of L of the two growths seems to be 
quite subtle and largely not yet understood at the mathematical level. The only rig- 
orous results in this direction are those obtained in [6] where, apart from logarithmic 
corrections, the appropriate lower bounds on T re i ax and T m ; x have been established by 
means of quite subtle test functions combined with the whole machinery of the Wulff 
construction. 

As far as upper bounds are concerned, they proved to be quite hard to obtain and the 
available results are still quite poor. In the case of homogeneous boundary conditions 
it was first shown in [16j that, for j3 large enough and any e > 0, 

^relax < exp ( cL 1/2 ~ 



for a suitable constant c depending on e and (3. Later such a bound was improved to 
exp(c^L log L) in [12]. When the inverse temperature f3 is just above the critical value, 
the only available result is much weaker (see [8] ) and of the form 

lim — log(T re i ax ) = 0. 

L— >oo Li 

Finally when /(c) = o~q the above bounds combined with some simple monotonicity 
arguments prove that, for any a > 0, 

Var+ {e tc f) < c/t a 

(where Var+ denotes the variance w.r.t. the plus phase 7r+ ) while the expected behavior 

is 0(e _v ^), see [TO]. 

We are now in a position to state our main results. 

Theorem 1.6. Let (3 be large enough and let L belong to the sequence {2 n — l} n£ ^. 
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(1) If the boundary conditions (b.c.) r are sampled from a law P which either 
stochastically dominates the pure phase 7r+ or is stochastically dominated by 



tt^ (see Section 2.2), there exists c = c{(3,e) (independent ofP) such that 

E||/i± -7rl<exp(-cL e2 / 16 ), (1.11) 

where ti = exp(cL £ ). In particular, 

P (T mix > t L ) < exp (-cL £2 / w ) . (1.12) 



(2) The estimates ( |l.ll[ )-( |i.l2 ) hold also ifP is stochastically dominated by ir c 



on 

one side of Al, and stochastically dominates 7r+ on the union of the other three 
sides. Similarly if the role of + and — is reversed. 

The most natural consequence of the above result is 

Corollary 1.7. Let (3 be large enough and let L belong to the sequence {2 n — l} n ^. 
Consider the square Al with b.c. r = +. For every e > there exists c = c{(3,e) < oo 
such that 

T mix <e cL \ (1.13) 

The same bound holds if the boundary conditions are + on three sides and — on the 
remaining one. Similarly if + is replaced by —. 

Remark 1.8. 



(i) In the proof of Theorem 1.6 and of Corollary 1.10 below, we need at some point 
some key equilibrium estimates which are proved in the appendix via standard cluster 
expansion techniques for values of (3 large enough. However, we expect those bounds to 
hold for every (3 > (3 C . Since this is the only part of the proof where the value of (3 comes 



into play, we expect Theorem 1.6 and Corollary 1.10 to hold for any (3 > (3 C . Let us 
also point out that, while we restrict for simplicity to the nearest-neighbor Ising model, 
we believe that our techniques can be generalized without conceptual difficulties to fer- 
romagnetic Ising models with finite-range interactions. In particular, cluster expansion 
results for large (3 are known to hold also in this more general situation. 

(ii) The restriction that L belongs to the sequence {2 n — l} ng pj is purely technical 
and it is a consequence of the iterative procedure we use. It would not be difficult to 



eliminate this restriction by somewhat modifying our iteration below (see Remark 3.12 



at the end of the proof of Theorem 3.2), but we have decided not to do this, in order to 
keep the presentation as simple as possible. 

(iii) The above results have been stated for the heat-bath dynamics but they actually 
apply to any other single site Glauber dynamics (e.g. the Metropolis chain) with jump 
rates uniformly positive ( e.g. greater than 5 > 0) as can be seen via standard comparison 
techniques [T7] . More precisely, if T m [ x and T ro i ax denote the mixing and relaxation 
times of the new chain, then there exist constants c, d depending on 5, (3 such that 
Tmix < c|A|T re i ax < c'|A|T re i ax < c'|A|T m i x ; the results we are after then follow since |A| 
represents a polynomial correction which is irrelevant in our case. 



(iv) Notice that in some sense our result (1.12 ) is not so far from optimality. Indeed, 



consider the distribution P such that r = + except for the boundary sites which are at 
distance at most L e from one of the corners of the box, where r is sampled from 7r+ . 
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Clearly P stochastically dominates 7r+ . Then, with P -probability exp(— cL £ ), r = — 
around the corners and, thanks to the results o/[T] ; T m i x > exp(cL £ ). 

1.4. Applications. It is intuitive that if the b.c. are all + (all — ) and we start from 
the all + (all — ) configuration, equilibration will be much quicker. Indeed, we have the 
following 

Corollary 1.9. Let (3 be large enough and r = +. For every e > there exists 
c = c(/3, e) > such that 

lim ||yu+ -7T-II =0, (1.14) 

L^oo 

where t\ : = exp(c(log L) £ ) . By a global spin flip the same results holds if + is replaced 
by -. 

Finally, here is the result about the decay of time auto-correlations for the infinite- 
volume dynamics in a pure phase: 

Corollary 1.10. Let (3 be large, let f(o~) = o~o and let p(t) = Var+ (e* £ /) be the time 
auto- correlation of the spin at the origin in the plus phase 7r+ . Then for any e > 
there exists a constant c = c{(3, e) such that 

p(t)<ce~^°^ 1/E . (1.15) 

2. Auxiliary definitions and results 

In this section we collect some more detailed notation that will be needed during the 
proof of the main results, together with certain additional auxiliary results that will 
play a key role in our analysis. 











L 




L 




> 

L 


A 


< * 

L 



Figure 1. The rectangle A and its enlargement El{N) 

2.1. Geometrical definitions. The boundary of a finite subset A C Z 2 , in the sequel 
denoted by dA, consists of those sites in Z 2 \ A at unit distance from A. Given a 
rectangle A C Z 2 and L 6 N, we denote by El(A) the enlarged rectangle obtained 
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from A by shifting by L units the Northern boundary upwards, the Eastern boundary 
eastward and the Western boundary westward (see Figure [T]) . 
Given e > (to be thought of as very small) and L G N we let 

R e L = G 1? : 1 < * < L, 1 < j < 

Similarly we define the rectangle Q L , the only difference being that the vertical sides 
contain now \(2L + l)2 +e ] sites. 

Notation warning. In the sequel we will often remove the superscript e from our 
notation of the various rectangles involved since it is a (small) parameter that we 
imagine given once and for all. 

2.2. Boundary conditions. A boundary condition r for a given domain (typically, a 
rectangle) is an assignment of values ±1 to each spin on the boundary of the domain 
under consideration. 

Definition 2.1. A distribution P of b.c. for a rectangle R (which will be Rl, Ql or 
a rectangle obtained by translating one of them by a vector v G 1?) is said to belong to 
T>{R) if its marginal on the union of North, East and West borders of R is stochastically 
dominated by (the marginal of) the minus phase ir^ of the infinite system, while the 
marginal on the South border of R dominates the (marginal of the) infinite plus phase 

7T+ 
"oo • 

The most natural example is to take P concentrated on the boundary conditions r 
given by r = — on the North, East and West borders, and r = + on the South border. 
In that case we will sometimes write tt r ' ,+ ' for the equilibrium measure in R, where 
we agree to order the sides of the border clockwise starting from the Northern one. 

2.3. The inductive statements. Here we define two inductive statements that will 
be proved later by a "halving the scale" technique. 

Definition 2.2. For any given LgN, <5>0,i>0 consider the system in Rl, with 
boundary condition r chosen from some distribution P. We say that A(L,t,5) holds if 

E\\fif -7T T \\<6 (2.1) 

for every P G D(Rl). 

The statement B(L, t, 5) is defined similarly, the only difference being that the rec- 
tangle Rl is replaced by Ql (in particular, P is required to belong to T>(Ql))- 

2.4. Censoring inequalities. In this section, we consider the Glauber dynamics in a 
generic finite domain A cZ 2 , not necessarily a rectangle. The boundary conditions r 
are not specified, because the results are independent of it. 

A fundamental role in our work is played by the censoring inequality proved recently 
by Y. Peres and P. Winkler: this says, roughly speaking, that removing (deterministi- 
cally) some updates from the dynamics can only slow down equilibration, if the initial 
configuration is the maximal (or minimal) one. 

First of all we need a simple but useful lemma: 



10 



FABIO MARTINELLI AND FABIO LUCIO TONINELLI 



Lemma 2.3. |22| Lemma 16.7] Let tt,h,v be laws on a finite, partially ordered proba- 
bility space. If v < n and v/ir is increasing, i.e. 

via) v(ri) 

-H>-TT 2 - 2 
7T{a) TT{r]j 

whenever a > r), then 

\\ v — tt|| < Ha* - tt||- (2-3) 

The result of Peres- Winkler can be stated as follows: 

Theorem 2.4. |22| Theorem 16.5] Let m E N, v := (ui, . . . , v m ) a sequence of sites in 
A, and let v' be a sub-sequence of v. Let ij,q be a law on such that ho/tt is increasing. 
Denote by fj, v the law obtained starting from fiQ and performing heat-bath updates at the 
ordered sequence of sites v. Similarly for /v . Then, 

\\fJ>v — Ti" 1 1 < ||/V ~ ^ll (2-4) 
and fj, v ^ jx v i. Moreover, [i v /ir and ^l v '/tt are increasing. 



It is easy to see that, if hq/tt is instead decreasing, (2.4) still holds, while the other 
statements become /v X fj,„ and [Av/n, /U^'/ 71 " decreasing. 

Here, "performing a heat-bath update at a given site v S A" simply means freez- 
ing the configuration outside v and extracting a v from the equilibrium distribution 
conditioned on the configuration outside v. 



Theorem 2.4 is proved in [52] in the particular case where is the measure con- 
centrated at the all + configuration, but the proof of the above generalized statement 
is essentially identical. Let us emphasize that such result is not specific of the Ising 
model but requires in an essential way monotonicity of the dynamics. 



From Lemma 2.3 and Theorem 2.4 we easily extract the continuous-time censoring 



inequality we need: 

Theorem 2.5. Let n £ N, = t < ti < . . .t n = T and Aj C A, % = 1, . . . , n. Let 

Ho be a law on S7a such that [jlo/tt is increasing. Let fx? be the law at time T of the 
continuous-time, heat-bath dynamics in A, started from (jlq at time zero. Also, let fi' T 
be the law at time T of the modified dynamics which again starts from no at time zero, 
and which is obtained from the above continuous time, heat-bath dynamics by keeping 
only the updates in Aj in the time interval [tj-ijtj) for i = 1, . . . ,n. Then, 

WfJ'T — "tt 1 1 < \\fpp — vr||, (2-5) 
and m -< fi' T ; moreover, — are both increasing. 



Needless to say, if instead ^q/tt is decreasing then all inequalities except (2.5) are 
reversed. 

Proof. Let m be the (random) number of Poisson clocks which ring during the time 
interval [0, T), and denote by Sj and Vi G A, % < m the times and sites where they ring. 
We order the times as Sj < Sj+i and of course Vi are IID and chosen uniformly in A. 
Define then w := ((v i, s\), . . . , (v m , s m )) and let [i w be obtained from //o performing 
single-site heat-bath updates at sites v\,V2, ■ ■ ■ ,v m (in this order) . Analogously, let w' 
be obtained by w by removing all pairs (vj, Sj) such that Vj ^ A^ where k is such that 
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Sj E [tk-i,tk), and fi w > be defined in the obvious way. For any realization of w one has 
from Theorem 2.4 that fi w ^ fi w i and that both fi w /"K and /J, w '/ir are increasing. Since 
fiT (respectively fi' T ) is just the average over w of fi w (resp. of fi w '), one obtains all the 
claims of the theorem (except (2.5)) by linearity. Inequality (2.5) comes simply from 



fJ'T ^ Mtj P^ us Lemma 2.3 and the fact that /xt/vt is increasing. □ 



We will need at various instances the following easy consequences of the above facts. 

Corollary 2.6. Let t > and assume that fi$ /it is increasing. Denote by fit the evolu- 
tion started from fi t =o = MO; o,nd by fif the one started from the maximal configuration 
+ . Then 

\\m-n\\ < \\i4-ir\\. (2.6) 



Proof. We know from Theorem 2.5 that fit/ft is increasing. Mo reov er, by monotonicity 

□ 



2.3 



of the dynamics fit ^ fif . The claim then follows from Lemma 

Corollary 2.7. Let j(t) = max (\\fif — 7r||, \\fij — 7r||). Then 

j(t + s) < 4rf(t)-y(s) W,s>0. 

Proof. Notice that \ \fif + s — n\\ = fi^ +s {A) — ir(A) where A = {a : fif +s (a) > vr(<r)}. 
Because of Theorem |2.5| the event A is increasing so that / := 1a — tt{A) is an increasing 
function (and of course 7r(/) = 0). Thus 

\W+s-*\\ = vt +8 { A )-< A ) 

= M t + 04(/)) 

= /4" (/£(/))-*(#(/)) 

< 2 7 (t)sup|^(/)| 

< 2 7 (t)max{| M +(/)|,|/i7(/)|} 

< 4 7 (t) 7 ( 5 ). 

Similarly for ff~ . □ 

2.5. Perturbation of the boundary conditions and mixing time. Consider a 
finite set A and two boundary conditions r, r. Let T m ; x and T m j x be the associated 
mixing times for the Glauber chain in A with b.c. r and f , respectively. Let M = 
max-tll^Hoo, Halloo}- 

Lemma 2.8. There exists a constant c independent of A,r, f such that 

< cM 3 |A|T mix . (2.7) 



Proof. Thanks to (1.10) and to the variational characterization of the relaxation time 
we get 

T mix < c|A|T relax < c|A|M 3 f relax < c|A|M 3 f mix 

where the third power of M comes from expressing the Dirichlet form, the variance and 
the local variances w.r.t. ir T in terms of those w.r.t. ir T . □ 
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Let now A C <9A, let ta be some configuration in Qa 5 let P be some distribution 
over the boundary conditions on dA and let P A be the distribution which assigns 
probability zero to b.c. r not identically equal to ta on A and whose marginal on 
dA \ A coincides with the same marginal of P. Notice that we can sample from P A by 
first sampling from P and then changing (if necessary) to ta the spins of r in A. If 
the pair so obtained is denoted by (t, f ) then the corresponding constant M satisfies 
M<M A :=e 8 ® A l 

Let d ± (t) = || fif — ir T \\ so that j(t) = max{d + (t),d (t)}. Similarly for d^(t),j(t). 
Lemma 2.9. With the above notation 

E( 7 (i))<e- M *+8E ( 7 (i)) 

where i = t/(c\A\ 2 M A ). 

Proof. Thanks to (|2TtJ) and fll~9j ), 

E ( 7 (i)) < e~ MA + P (T mix > t/M A ) < e~ M * + P (f mix > i/(c|A|Af£)) 
= e- M *+p(f mix >\A\i) 



Notice that, for any s > 0, T m ; x > s implies that there exists some starting configuration 
a for which the variation distance of its distribution at time s from the equilibrium 
measure ir T , call it d a (s), is at least l/(2e). However, using the global monotone 
coupling of the Glauber chain, 

d*(s) < P ( % + < f + Vs'') < E Mvt^(x) = +)- nVs''(x) = +)] (2-8) 



and therefore 



< |A| ld+( 8 ) + d~{s)j < 2\A\j(s) (2.9) 
T mix > |A|t) < P U(\A\i) > — L 



Thanks to Corollary 2.7 7 (t) < (4 7 (t )) L * /toJ so that 



P (m\t) > ^) < P > I) < 8E ( 7 (0) • 



□ 



Let us remark for later convenience that, exactly like in (2.8), one proves that 

sup |K -vr T || < 2|A| 7 (i). (2.10) 

a 

With the same notation the following will turn out to be quite useful: 

Corollary 2.10. Let R L = R\ and let P G V(R L ). Let also A C dR L be such that 
L 3e < |A| < 2L 3e . Assume that E A (\\fif - vr T ||) < 6 for every P G V(R L ). Then the 

statement A(L, t! , 5') holds true with 5' = 85+e~ e and t! = te for some constant 
c > independent of A and ta- Analogously A(L,t,S) implies E A (||/^7 — tt t ||) < S' . 
Similar statements hold if we replace Rl by Ql and A(L, t' , 5') by B(L, t! , 5'). 
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3. Recursion on scales: the heart of the proof 

This section represents the key of our results. We will inductively prove over the 
sequence of length scales L n = 2 n+1 — 1 that the statement A(L n , t n ,5 n ) and its analog 
B(L n ,t n ,6 n ) hold true for suitable t n ,S n (see Theorem 3.2 below). In all this section 
e > is fixed very small once and for all. Accordingly, for any L E N, Rl = -Rf, and 
similarly for Ql- Finally c, d will denote positive numerical constants whose value may 
change from line to line. 

First we give a rough estimate which provides the starting point of the recursion: 

Proposition 3.1. For every (3 there exists c = c{(5) such that for every L E N the 
statements A(L,t,e~ te ° L ) and B(L,t, e~ te cL ) hold. 



Proof. From rough estimates on the spectral gap |16| Corollary 2.1] and (1.10), one has 
that 

T mix < e cL (3.1) 

uniformly in the boundary conditions r and in L £ N, both for Rl and for Ql. Applying 
(1.9) with e = l/(2e), the claim is proved. □ 

Theorem 3.2. For every [3 there exist constants c,c' such that: 

(1) if A(L,t,5) holds, then also B{L,2t,5\) does, with 

Si =5 1 (L,5,t) =c{5 + e- c ' L2e + L 2 e' c ' l °^) . 

(2) If B(L,t,5) holds, then also A(2L + 1, t2, 62) holds, with 

t 2 = t 2 (L,t) = e cL3£ t (3.2) 

and 

5 2 = 5 2 {L,8) = c{5 + e- c ' L ' ie ). (3.3) 

Assuming the theorem we deduce the 

Corollary 3.3. There exist c,c' > such that the following holds. For every L G 
{2 n — l} n£ N there exists 

A{L) < exp (-c'i £2 ) (3.4) 
such that A (L, t, A(L)) holds for every t > e cL,ie . 

Proof. Note that if one iterates j times the map x \—* 2x + 1 starting from x = 1 one 
obtains 2 J+1 — 1 =: Lj. Assume now that L = L n for some large n and set no := L en J > 
so that (l/c)L £ < L no < cL £ . 

From Theorem 3.2 one sees that it is possible to choose c, c' > such that 

A(Lj,tj,5j) => A{L j+ i,t j+ i, 5j + x) (3.5) 

with 



and 



t j+ i = 2 tj e cL f (3.6) 



8 j+1 = c (Sj + e~ c ' L f + L) e - c ' l °^) . (3.7) 
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A 


\{2L + l) 1 ^+ e ] - ^ 1/S+E ] 




' 2\L$ +e ~\ - [(2L+l)i +g ] 

r 


B 



Figure 2. Ql and its covering with the rectangles A, B 



Let 



cL 3£ 



so that, thanks to Proposition 3.1, A(L no , t no , 5 no ) holds with 



J no 



exp 



cL 3£ 



Then, applying (3.5) n — uq times, one obtains the claim A(L,T(L),A(L)) with 

T(L) : = 2 n - n ° e cE ?=™o L f < e ci3e 

and 



A(L) < L c 5{n ) + ( e- c ' L "° + e 



-c'log(t„ ) 



< e 



-cL e 



for a suitable constant c, where we used the rough bound (cf. (|3.7|)) 



(3.8) 

(3.9) 
(3.10) 

(3.11) 
□ 

E([| M +-tt1) <<5l (3.12) 

Observe that can be seen as the union of two overlapping rectangles A and -B, 
where B is just the basic rectangle Rl and A is obtained by shifting B to the North 
by \{2L + 1)V2+*1 - [L 1 / 2+£ ] (see Figure §. 

Let now /4 denote the distribution at time 2t of the dynamics started from the all 
+ configuration and subject to the following "massage": in the time interval [0, t) we 
keep only the updates in A, at time t we increase all the spins in B to +1 and in the 
interval (t, 2t] we keep only the updates in B. 



The statement for every t > T(L) then follows from Corollary 2.7 



3.1. Proof of Theorem [jTg} part (1). 

i) We begin by proving that for every distribution P 6 T>{Ql) one has 



Lemma 3.4. 



|/4 -^11 < II Aai 
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Proof. Let fi^t denote the distribution at time It of the dynamics started from the all 
+ configuration and subject to the following "censoring": in the time interval [0,t) 
we keep only the updates in A and in the interval [t, 2t] only the updates in B. By 

^pr is increasing. Moreover fi^ di fitt wrnc h combined with Lemma 



2.5 



2.3 



Theorem 

proves the result. □ 

In order to better organize the notation we need the following: 
Definition 3.5. We let 

(a) v\ be the distribution obtained at time t after the first half of the "massage". Clearly 
v\ assigns zero probability to configurations that are not identical to + in A c ; 

(b) ^2 be the distribution obtained from the second half of the censoring starting (at 
time t) from a configuration equal to + in B and to a in B c . Clearly v% assigns 
zero probability to configurations that are not identical to a in B c ; 

(c) vr^ + := tt t (- | a a- = +); 

(d) Tr T jf := vr T (- 1 a B c = rj); 

(e) tt t ~ (resp. tt t,+ ) be the Gibbs measure in Ql with minus (resp. plus) b.c. on its 
South boundary and r on the North, East and West borders. 

With these notations the distribution jj^t ls written as 
Notice that also the Gibbs measure ir T has a similar expression, namely, 



Therefore 



where 
Clearly 



V 

V V 

= \\u\ - 7r^' + || B c + ||7 - -71- 1| (3.13) 
7(77) := 7t^ + (t/ Bc ) vl BC {r]). 



It - A < ^~{\W V 2 BC - *W BC II) + hT - * T \W + h T - ^Wbc 



In conclusion 

E(||/4-7rl) <E(|h-7r5 + [|flc 
+E (tt^- (||^ c - tt^ II)) + E (||tt^ + - Trifle) + E (K - ■ (3- 1 I ) 



By assumption the first term in the r.h.s. of (3.14) is smaller than 5. Next we analyze 
the second term. In this case, if we denote the four boundary conditions around B, 
ordered clockwise starting from the North one, by 7-1,72,73,74, then their distribution 
P is given by 

P~(Tl,T2,73,T4) = P(t2,T3,T 4 )E (7T T ~(ti) | r 2 ,r 4 ) . 
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Notice that the marginal of P~ on T3 coincides with that of P and therefore stochasti- 
cally dominates the corresponding marginal of 7r+ . It remains to examine the marginal 
on (ri,T2,T4). Let / be a decreasing function of these variables and observe that, as 
a function of the boundary conditions on the North, East and West sides of Ql, the 
average TT T ~(f) is also decreasing. Therefore, since P £ D(Ql), 

E-(/) = E > 7T- > 7T- (/) (3.15) 

i.e. P" G £>(£). Therefore 

e (^-(ii^r - ^ bc \\)) = e- oi*r - ^ s d < & 



The third and the fourth term in (3.14) can be bounded from above by essentially the 
same argument which we now present only for the fourth term. Clearly, for any choice 
of the boundary conditions r, 7r r,_ < ir T . Therefore 

E (K " 7T T '-\\ B c) <Y, E = +) - = + )) • 

x<=B c 

Claim 3.6. There exists c = c(/3, e) > such that 

E (n T {a x = +)- n T '-(a x = +)) < e~ cL * e (3.16) 

for every x G B c . 

Proof. Let T denote the event that in B there is a *-connected chain {i.e. either the 
Euclidean distance between two consecutive vertices v,v' of the chain equals 1, or it 
equals y/2 and in that case the segment vv' forms an angle ir/4 with the horizontal 
axis) of — spins which connects the East and West sides of B. By monotonicity, 



and therefore 



n T (a x = + I T) < vr T <->, = +) (3.17) 
tt t (o- x = +)- 7r T '-(a x = +)< tt t (T c ). 

By monotonicity 

E = +) - ^~{a x = +)) < Evr-(r c ) < tt" (^ r ' + (r c )) 

where we recall that the superscript + means that on the South border of Ql the b.c. 
are all plus. Let tt^o ' ' be the the minus phase measure -n - ^ conditioned to have all 
minuses on the North, East and West borders of the enlarged rectangle El{Ql) (see 
Figure [3j . Standard bounds on the exponential decay of correlations in the minus phase 
(see for instance [20] or (2U Chapter V.8]) prove that 

^ (vr r ' + (r c )) < 7&'-> K' + (r c )) + e~ cL (3.18) 

for some constant c > 0. If we now add extra plus b.c. on the whole horizontal line 
containing the South boundary of Ql and denote by itio ' the corresponding Gibbs 
measure then, by monotonicity and DLR equations, we obtain 

t&'-> K' + (r c )) < t4<-' +) K< + (r c )) = ,&.-.+) (r c ) . (3.19) 

Notice that irio ' is nothing but the Gibbs measure ^e'^Ql) m ^ e rec t an gl e El{Ql) 
of Figure [3j with + b.c on the South border and — b.c. on the rest of the boundary. 
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Figure 3. The rectangle Ql (thick line) and its enlargement El{Ql) 
(narrow line), with the b.c. of 7rL ' 



Next, note that the event T c implies that the unique open Peierls contour 7 (see 
definition in Appendix |Al crosses the horizontal line containing the South border of A, 



and we will prove in Appendix 1.2 that 



7T 



El{Ql) 



(7 reaches the height of the South border of A) < e 



cL 2e 



(3.20) 



The intuition for (3.20 ) is that the open contour 7 behaves like a one-dimensional simple 
random walk starting at the origin and conditioned to stay positive and to return at 
time L to the origin: the probability that before this time it goes at distance of order 
L l l 2+e from the origin is smaller than exp(— cL 2e ). □ 



Altogether we have obtained 

E||/4 -7r r || <2J + e" 



ii) Now we consider the dynamics started from the all — configuration and we prove 

E||^2t - 7r T || < ft. (3.21) 

By Theorem 



2.5 



\fi2t ~ 7r ' r || — 1 1 Asm — ^ll where this time Jl^t denotes the distribution 
at time It obtained by starting the Glauber dynamics from the minus initial condition 
and performing the following "massage" (the reverse of the previous one) : in the time 
interval [0, t) we keep only the updates in B, at time t we reset to — all the spins in 
A and in the time interval [t, 2t] we keep only the updates in A. In order to keep the 
notation as close as possible to that of the previous case where the starting configuration 



was all pluses we redefine 

Definition 3.7. 

(a) tt t b 

T.r 



(b) 



= 7T T (. 
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(c) v\ is the distribution obtained after time t and is that obtained in the second 
time lag t starting from the configuration equal to — in A and to a in A c . 



With these notations the same computation leading to (|3.14|) gives 
E||/2m 



TT 



TT 



B 



\ A c+E7r T (\wr c 



TT 



7^||)+E||7r T R 



B 



" 7T ||^c 
cL 2e 



(3.22) 



The first and third in the r.h.s of (3.22) are smaller than 5 and e _cij "" respectively 
by essentially the same arguments as before. It remains to analyze the second term. 
Notice that 



TT 



MA? 



TT 



<5> r [*2 A >. = W Ae (^ = -))] 



(cr x 



tt t {o x 



Given x £ A and I G N, let Kg be the intersection of A with a square of side 2£ + 1, 
centered at x. Monotonicity implies that 



(cr x =-)< V%f{<T x = ~) 



(3.23) 



where v2ff denotes the distribution at time t obtained by the dynamics in Kg, started 
from all — , and with b.c. which are all — except on dKi n dA where the b.c. remain 
either r (on the North, East and West border of A) or r]A c (on the South border of A). 
Let 7tT'' Mc be the equilibrium measure of this restricted dynamics. Then, 



< 



.VAC 
^21 



{<T X = ~) 



TT 



MA? 



t,Vac 



(o- x ■ 
(o- x 



- 7T T (a x = -) 
= -) - n T (a x 



iT T (a x 



-)] 



where in the last inequality we used (3.1). If we now average first with respect to tt 7 
and then with respect to P we claim that 



Claim 3.8. On has for some c > 

E(^(^> 2 



-))-*•>« = -)) 



E (tt t [tt™ c (a, = -) - tt t /> ac {cr x = -)])< e 



-ci 



(3.24) 
(3.25) 

and the left-hand 

side of fl3T24~l ) equals 0). 

Assuming the claim it is now sufficient to choose £ = [~(l/c)(logi — log log t)] to 
conclude that 



(It is clear that if £ is so large that -fQ = A, then 7rI' T?AC 



TT 



t,Va c 



TT 



t,Vac I 



< L 2 e- C ' l °^ 



for some d > 0. 



(3.26) 
□ 



Proof of Claim 3.8 Let Y be the event that x is separated from dK^ D ^4 by a *- 
connected chain of minus spins. By monotonicity, for any t]a c i 



TT 



T,VA C 



(cr x 



T) > TT 



T,r]A c 



(cr x = -) 
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and therefore it is enough to show that 



E(^7-(r c )))=E(^(r c ))< e - 



-ci 



The rest of the proof is now very similar to that of Claim 3.6 Apart from an error e 



-cL 



we can replace E (7r r (T c )) by k e Wq \ (r c ), where ^e l (qI) ls Gibbs measure on the 
enlargement El(Ql) (see again Figure [3]ab ove) with plus b.c. on the South border and 
minus b.c elsewhere. In turn, thanks to the fact that the event T c depends only on the 
spins in A, we can replace ^e'^Ql) ^Y the Gibbs measure 7r E L (Q L ) on ^ ne same region 
but with homogeneous minus b.c. by paying an error smaller than e~ cIj2e . Finally, 
again by monotonicity and standard correlations decay bounds in the pure phase, 

.(r c )<7r-(r c )< e - rf 



7T 



El(Ql) 



for some c > 0. 



□ 



3.2. Proof of Theorem 3.2 part (2). Thanks to Corollary 2.10 and apart from 
the harmless rescaling t i— > tf = e cLiE t and S \—* 5' = d8 + e~ cL e for some constants 
c, c' > 0, we can safely replace the distribution P over the bo unda ry conditions outside 
R2L+1 with the modified distribution P A (defined in Section 2.5), where A = {(i,0) G 
dB.2L+il I* — L\ < L 3e } and the pinned configuration ta is identically equal to —1. In 
other words it is enough to prove that E A (\\fJ^ t , — vr T ||) < cS'. 



i) As before we begin with the case where the dynamics in R2L+1 is started from all 
pluses. Let now (see Figure [4]) 

A = Q L + ([L/2\,0) 

B = {Q L }U{Q L + (L + 1,0)} 

C = eR 2L +i; i = L + l}. 



so that i?2L 



+1 



BUC and B(lC 



B 



B 



Figure 4. R2L+1 and its covering with A, B, C. In bold the set A 



By Theorem 



2.5 



7T 



< 



\P-tt' 



7T T || where, as before, the tilde indicates that 



the following "massage" has been applied: in the time interval [0,tf) we keep only the 
updates in A, at time t 1 we increase to +1 all the spins in B and in the interval (t', 2t'\ 
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we keep only the updates in B. Notice that the dynamics in B in the time lag (t' , 2t'\ 
is a just a product dynamics in the two copies of Ql, in the sequel denoted by B\ and 
-B2, whose union is B, with boundary conditions r on dBndR2L+i and some boundary 
conditions on C generated by the dynamics in A in the first time lag [0,t']. 

Definition 3.9. We define 

(a) v\ as the distribution obtained at time t' after the first half of the censoring; 

(b) v*2 as the distribution obtained from the second half of the censoring starting (at 
time if) from a configuration equal to a in C and to + in B. Clearly Uq, assigns 
zero probability to configurations that are not identical to a in C; 

(c) tt^ + := 7r T (- 1 ua c = +) and similarly with + replaced by —; 

(d) ti t ^ c := ir T (-\a c = rj c ); 

(e) 7r r '~ (resp. ir T ' + ) as the Gibbs measure in R2L+1 with minus (resp. plus) b.c. on 
its South boundary and r on the North, East and West borders. 



By proceeding exactly as in the proof of statement (1) we get 

- ""I < ft -7T T || 



< \\vi-KT\\c + K T >-{\wf 



B 



+ \\tTa — it t \\c + ||vr T ' — vr T ||c 



(3.27) 
(3.28) 



and 



E A (||M+,-vrl) <E A (||^-^ + || c 



-E A (7r T '- 



7T 



7T 



7T ||c 



+ E Z 



tt t '-\\c). (3.29) 



By assumption and thanks to Corollary 2.10 if we perform a global spin flip we see 



that the first term in the r.h.s. of (3.29) is smaller than 5'. As far as the second term 
is concerned we observe that the distribution P A >~ of the boundary conditions (r, rjc) 



given by P A - (t, rj c ] 
P-{t, VC ) = P(rK 
belongs to V(Bi), i ■ 
||// (8) v — // (g) < I 



P (t)tt t ' {rjc) coincides with the A-modification (P ) of 



(rjc)- The same argument as in (3.15) shows that the latter 



1,2, so that (via Corollary 2.10 and the immediate inequality 
,/ 1 



fi - /i || + \\u 



the second term is smaller than 25' . 



We now turn to the more delicate third and fourth term in the r.h.s. of (3.29). Since 



they can be treated essentially in the same way we discuss only the third one. As usual 
we write 



E- 



7T 



7T ||C 



< Y, eA K'Vz = +) " ^ r (^ = +) 



(3.30) 



x&C 



Let r be the event that in A there exist two *-connected chains of minus spins, one to 
the left and the other to the right of C, connecting the South side of A to its North 
side. By monotonicity 



7T 



T,+ 



(a x = + I r) - w T (a x = +) < 



so that 

vr^ + (a x = +) - 7r T (a x = +) < vr^ + (r c ). (3.31) 

Let now A = 1 < i < L, 1 < j < 2\{2L + l) 1 / 2 ^]} so that A consists of just 

two copies of A stacked one on top of the other. Then, using monotonicity together 
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with the standard exponential decay of correlations in the minus phase tt^ (see e.g. 
Q3.18D ) we get 

E A (irTtn) < e~ cLl/2+£ + ^' +A \n (3-32) 



A \ J J — ^ ' n a 

where the superscript (— , +, A) indicates the b.c. which is — on the union of the North 
boundary and A, and + on the rest of dA. The key equilibrium bound we need at this 
stage is the following: 

Claim 3.10. There exists c> such that vr^"' + ' A) (r c ) < e~ clJiE . 



Putting together the bounds we got on the various terms in (3.29), we have proved 
E — vr T || < c5' as wished. 

The proof of the claim is deferred to the appendix but intuitively the argument goes 
as follows. Under the boundary conditions ( — ,+,A), for any configuration a £ 0,^ 
there exist exactly two open Peierls contours 71,72 with two possible scenarios: 

(a) 71 joins the two upper corners of A and 72 the two ends of the interval A; 

(b) 71 joins the left upper corner of A with the left boundary of A and similarly for 72. 



If we recall the definition of the surface tension (1.2), the ratio between the probabilities 
of the two cases is roughly of the form: 

e -/3T /3 (ei)(L+2L 3£ )+2/3r (3 (e)D 

where D is the Euclidean distance between the left upper corner of A and the left 
boundary of A and 9 is the angle formed by the straight line going through these two 
points with the horizontal axis. Clearly 9 ra 0(L~2+ £ ) and D m L/2 - L 3e + 0(L 2e ). 
Therefore case (b) is much more likely than case (a). 

Remark 3.11. Notice that it is exactly the presence of the positive correction 0(L 2e ) 
in D that forced us to take the length of A to be L 3e . 

Once we are in scenario (b) the most likely situation is that neither 71 nor 72 touch C 
(otherwise they would have an excess length of order L 3£ ) and the desired bound follows 
by standard properties of the Ising model with homogeneous boundary conditions. 

ii) The proof of E A ||/^ t , — ^ T \\ < cS' is identical, modulo the obvious changes, 
provided that we redefine the "massage" of fxZ/ as the censoring in A, B plus the 
resetting at time t' of the spins inside B to the value —1. A minor observation is that 

in this case, for the smallness of the term E A — vr^' - , we do not need anymore 
the global spin flip that was necessary for the dynamics started from all pluses. □ 

Remark 3.12. As we said at the beginning, in order to keep the focus on the main 



ideas of the method, Theorem 3.2 has been given in the restricted setting in which the 
length scales are of the form L n = 2 n — 1. However it should be clear by now that 
the case of arbitrary length scales can be dealt with in a very similar way. A possible 
solution requires a slight modification of the definition of the two inductive statements 
A(L,t,S),B(L,t,S). 

Let Tl (respectively Ql) be the class of rectangles which, modulo translations, have 
horizontal base L and height H S [L? +e , (2L) 2 +e ] (resp. horizontal base L and height 
H £ [(2L) 2+ e ; (4L) 2 +e ]J. Notice that any rectangle in Ql can be written as the union 
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of two overlapping rectangles in Tl such that the width of their intersection is still 
0(L 1//2+£ ) (as in Figure^). Moreover for any n large enough and any L £ [L n+ \, L n+ 2) 
there exists L' £ [L n , L n+ \) such that any rectangle A in Tl can be written as the union 
of three sets A, B,C (as in Figure^! where A £ Qy, B consists of two disjoint rectangles 
in Qy and C = A \ B satisfies dist(C, A c ) = 0{L) and has horizontal width O(l). 



We then say that A'(L,t,5) (B'(L,t,5)) holds if (2.1) is valid for every rectangle in 

continues 



J~L (i n Ql)- It is almost immediate to check that part (1) of Theorem 3.2 



to hold with this new definition. Part (2) can be modified as follows. If B'(L' ,t,5) 
holds for every L' £ [L n , L n+ i) then A'(L, t2, 62) holds for every L £ [L n+ i, L n +2) with 



t2 = e c23 " £ t and 62 = c(5 + e 
same as that given above. 



The proof of the new version is essentially the 



4. Proof of the main results 



In what follows we will prove Theorem 1.6 and Corollaries 1.9 and 1.10 Notice that, 
for any A C Z 2 , any boundary conditions r and any starting configuration a, — tt t \\ 
is invariant under the global spin flip r — r and a 1— > —a. Therefore it will be enough 
to prove only "half of the statements" . 



4.1. Proof of Theorem 11.61 Recall that 

ti := exp(cL £ ) 

for some chosen e > small, and let e' := e/4. We assume throughout this section that 
L £ {2 n - l} n£ N- 

4.1.1. Mixing time with "approximately (—,—,+,—)" boundary conditions. First we 
prove ( 1.11 )-( Jl.l2 ) when the b.c. r is sampled from a law P which is dominated by 71"^ 
on the union of three sides of A^ and dominates 7r^ on the remaining side (e.g. the 
South border). 



One sees from (2.10), the definition (1.8) of mixing time and the Markov inequality 
that ( 1.11| ) implies (1.12), so we are left with the task of proving ( |1.11[ ). This is an 
almost straightforward generalization of the proof of point (1) of Theorem 3.2 and 
therefore some steps will be only sketched. 

For definiteness, we assume that the LxL square Ax, we are considering is {(xi, X2) £ 
Z 2 : 1< xi,X2 < L}. Consider first the evolution started from the + configuration. 
For i > let 



hi 



L V2- 



+ i 



(2L+1 



1/2+e' 



L 1/2H 



To avoid inessential complications, assume that there exists k £ N such that h^- 
Of course, 

L l/2-e' 



21/2+e' _ 1 ' 



(4.1) 

= L. 
(4.2) 



Let A l L be the rectangle of height hi whose base coincides with that of A^, so that in 
particular A^ _1 = A^. We will prove by induction at the end of the present section 
that 
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Lemma 4.1. The following holds for i = 0, ...,k — 1. Let i/ie b.c. r around the 
rectangle A l L be sampled from a law P which dominates 7r+ on the South border and is 
dominated by 7r~ on the union of West, East and North borders. Then, 



E\\n 



7T» 



<(l + ^-^' )2 =(l + z)e-^ 2/16 , 



(4.3) 



"(i+l)t L /k 

where fj^' 1 is the evolution in A' l L started from +, 7rTj is its invariant measure and c 
depends only on (3 and e. 

If the Lemma holds, it is sufficient to apply it for i 
exp(— aZ/ 2 / 16 ) as wished. 

It remains to show that 



k— 1 to see that E||/^ — 7r r || < 



El 



7T 



< e 



By Theorem 2.5 and (the analog of) Lemma 



3.4 



-A < Win, 



(4.4) 

•7r||, where this time 



fi t is the dynamics in obtained via the following "massage": in the time interval 
[0,£l/2) we keep updates only in B := = {(xi,X2) G Aj, : x 2 < [L 1 / 2 ^']}, at time 
t L /2 we set to - all spins in A := {{xi,x 2 ) G A L : x 2 > {(lfflL 1 ^]} and in (t L /2,t L ] 
we keep updates only in A. In analogy with Definition |3.7| we introduce the 

Definition 4.2. We let 



(a) 7T 



B 



7T T (-|<7 B c 



= -); 
v); 



(b) it t J? :=it T {-\o A c 

(c) v\ be the distribution obtained at time t^/2; 

(d) v% be the distribution obtained at time tj_,, starting at time t^/2 from a in A c and 
from — in A. 



Then, in analogy with (|3.22|) one finds 
E||/V, 



7T 



7T, 



U c +Evr T 



,VAC 



7T 



+ Ellvr 



B 



^r T II 



(4.5) 



From Corollary 



3.3 



one sees that the first term is smaller than exp(— cL £ Z 16 ) (note 
that £l/2 3> exp(cL 3s )). The last term in (4.5) can be bounded by exp(— c'L 2e ) (the 
proof is essentially identical to the proof of the upper bound on the last term in ( |3.22 ) ) . 



Finally, proceeding like for the second term in (3.22), one sees that 

E 7T T 



u VAo _ r.^lh < L 2 e -c'log(t i /2) + -<M* <<; e -c'L* 



Altogether, we proved (4.4) and the proof of (1.13) is complete 



(4.6) 
□ 



Proof of Lemma 



4.1 



is just Corollary 3.3 (note that A° L 



Let for simplicity of notation tt t := 7rTj . For i = the claim 



Assume that the claim holds for i 



We define the following three disjoint rectangles (see Figure 5b: A := A l L \ A 



1. 

C 

is the rectangle whose South border coincides with t hat of A^ and whose height is 
^l/2+e'-^ and B ._ A i. \(A\j C). By Theorem 

< 



(h 



3.4 



m 



2.5 



and (the analog of) Lemma 
one has [|^ + x)tx,/fc ~~ 7I " T H — \\&7i+l)tz,/k ~ wnere the "massage" in /i^ consists 
seeping only the updates in yl U B in the time interval [0, t^/h) and in B U C in 
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h.-h. 



1/2+ E' 
h-L 



At 



Figure 5. The rectangle A l L and its decomposition into A,B,C. 

the time interval {ti/k, (i + Vftz/k], and setting to + all spins in B at time ti/k. In 
analogy with Definition |3.5| 



Definition 4.3. We let 

(a) v\ be the distribution obtained at time t^jk, which assigns zero probability to con- 
figurations which are not all + in C U B; 

(b) v% be the distribution at time (i + l)ti/k, starting at time ti/k from a in A and 
from + in B U C ; 

( c ) ^aub '■= 7r r (-kc = +); 

( d ) ^buc '■= n T (-\VA = rj); 

(e) tt t ~ be the Gibbs measure in A l L with — b.c. on its South border and r on the other 
borders. 



One has then 



and 



n T (rj) = -K T (r] A )ir T /^ A c (r]). 



(4.7) 
(4.8) 



In analogy with (3.13) 

H^(i+l)ii/fc 

where 



- 7T T || < \\l>i - W^jbWa + \\l ~ 7T T | 
7(77) := TTAuB^A^iv)- 



7T 



T,VA I 

Bud 



As a consequence, using (4.8), 

W^(i+l)t L /k - \Wi - ^aubWa + vr r '" 

+ll 7r J 4us ~~ w IU + If ~ ^ IU- 

Now we can take the expectation with respect to P. First of all, we have 



(4.9) 



(4.10) 



thanks to Corollary 3.3 because AUB is a translation of the rectangle R £ L which appears 



in the definition of the claim A(L, t,5). As for the P-expectation of the third and fourth 
terms, it is upper bounded by exp(— cL 2e ) (the proof is essentially identical to that 
of the upper bound for the third and fourth term in ( 3.14[ )). Altogether, the average 
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of the sum of the first, third and fourth terms is upper bounded by exp(— cL e2 / 16 ). 
Finally, in order to bound the P-expectation of the second term we need the inductive 
hypothesis. Indeed, we can say that 



E 7T T ' 



,VA 



< ie 



'2 -*BUC?llj ~~ 

(which concludes the induction step) if we prove that the marginal on the union of 
North, East and West borders of BUC of the measure E~ := E-7r r '~(-) is stochastically 
dominated by ir^. Indeed, if (ti,T2,T4) is a generic spin configuration of the North, 
East and West borders of B U C and / is a decreasing function, using monotonicity a 
couple of times one gets 

E-(/) = Ett t '-(/) > TP" (tt t '-(/)) > TP" (/), (4.12) 

which proves the desired stochastic domination. □ 



4.1.2. Mixing time with boundary conditions dominated by tr^. Here we prove ( |1.11[ ) 
(and therefore, via Markov inequality and (2.10), we obtain (1.12)), when the law P of 



r is dominated by tt^ (or, by spin-flip symmetry, when it dominates 7r+ ). 

We begin with the evolution starting from the + configuration and we recall that 
Ai = {l,...,L} 2 . One has by monotonicity tt t ^ fif, and therefore 



E||/4 



7T 



E E K 



+ 



ieA, 



+)). 



+)) (4-13) 



where A L := {(i,j) S : j < L/2} and A^ := A^ \ A L . We will show that the sum 
S\ is small, and 52 can be dealt with similarly. 



Recall that K l L and k were defined in Section 



4.1.1 



and observe that A^ 3 ^ 4 ^ is a 



rectangle whose base coincides with that of A^, and whose height is h ~ (3/4) L (cf. 



(4.1 )-(4.2)). Then, thanks to Theorem |2 . 5 1 (or actually by monotonicity), we know that 



/i? ■< fit , where fit is the censored dynamics in which only updates in A^ 
retained. One has therefore 

St < Yl E = +) - ^(^ = +)) 



L(3/4)fcJ 



arc 



(4.14) 



xeA, 



7T 



+ E||vr T 



7T 



r .+ l 



where tt t,+ is the invariant measure of jlf , i.e. 



7T 



7T 



+ 



Since the North border of A L is at distance approximately L/A from the North border 



of A^ 3//4 ^ fc -', the last term in (4.14) is easily seen to be upper bounded by exp(— c'L) 
(the proof of this fact is essentially identical to the proof of the upper bound for the 
last two terms in (3.14)). As for the first term, Lemma |4 . 1 1 ( applied with i = L(3/4)/cJ) 
shows that it is upper bounded by exp(— c'L £ Z 16 ). This is because the evolution jlf sees 
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b.c. + on the North border of A 



L(3/4)fcJ 
L 



, and t (sampled from P which is stochastically 



dominated by tt^) on the remaining three borders. Altogether, we have shown that 



7T 



< e 



_ C / L ^/16 



Next, we look at the evolution started from all — . Given a site x £ Al and £ £ N, 
let Ki be the intersection of Al with a square of side 2£+l centered at x. We let /J^~ t 
be the dynamics in Kg with — initial condition and with b.c. — except on dKg n 8Al, 
where the b.c. is r. The invariant measure of such dynamics is denoted by vr]^~. Since 
:< tt t , we have 



Elk" 



7T 



< 



< 



< 



E 

x£A L 

E [ E to 

x&Al 



-) - E7r T (a x = -)] 

= -) - Kk~((T X = 



+ e 



-cl 



E 



E 



7T 



A, 



+ e 



(4.15) 
(4.16) 
(4.17) 



The "error term" exp(— ctj comes from comparing 'EiTT T {a x = — ) and E7r^ {a a 



(see the proof of Claim 3.8 for very similar arguments). We know from |16[ Corollary 
2.1] that T^r K < e rf , uniformly in r. Therefore, from (1.9) and choosing t = tj, and 
I = [i(logt — log log i)] « L £ , one gets 



7T 



< e 



(4.18) 
□ 



4.2. Proof of Corollary 1.9 , We restart from ( |4.17 ), which in the case of r 



IT 



< \A L \e~ ce + ^ 



xeA L 



7T 



- gives 
(4.19) 



where it , /j, k t and ir K are just tt 1 ',//^^ and 7r^ respectively, in the specific case 
r = — . Now we use the extra information that the mixing time T~ ix K of the dynamics 



f^Kt t ^ S a ^ mos t exp(c'^ e ), as follows from (1.13). We choose i to be the smallest integer 
in the sequence {2™ — l} n gN such that cl > 31ogL, so that the first term in the r.h.s. 



of (4.19) is smaller than 1/L. Taking t\ := exp(c(log L) £ ), one has from (1.9) 



lA£> u ~ *k,\\ < e tl/T °**>Kt < exp[-exp(c(logL) £ - c'l e )} « 1/\A L \ 



(4.20) 



if one chooses c suitably larger than d (recall that we chose 
corollary is proved. 



O(logL)) and the 
□ 



4.3. Proof of Corollary 1.10[ This is rather standard, once (1.13) is known (cf. for 
instance Theorem 3.2 in [16J or Theorem 3.6 in [7]). Clearly, it is sufficient to prove 
the result with / redefined as /(c) := (co + 1) which has the advantage of being 
non-negative, increasing and with support {0}. Consider a square Je C Z 2 with side 
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2£ + 1 £ {2 n — l} n gN an d centered at 0. By the exponential decay of correlations in the 
pure phase 7r+ , 

\*W)-*%U)\<c<r dt - (4.2i) 

Moreover, by monotonicity, for every initial configuration a of the infinite system 



< 



(e tL f)(a)< (e**/) (a) 



(4.22) 



and the right 
of Section 



1.2 



rand side is an increasing function of a; in accord with the notations 
Cj t denotes the generator of the dynamics in Ji with + boundary 
conditions on dJi (its invariant measure is of course 7r^ ) and C is the generator of the 
infinite- volume dynamics. One has then (using once more monotonicity) 



7T ' 



[{e tC ff 



< IT 



Ji 



which, together with (4.21), gives 

p(t) = Vari {e tc f) < Va V 



-c'i 



By (1.6), one has that 



tc 



-2tgap- | ) 



Var w + (eT'tf) < Var w + (/)e 
with gap j the spectral gap of £j . From the inequality 



gap > 



1 



(cf. (1.10)) and (1.13), one deduces that for every e > 

-2te" rf£ 



Now letting 



Var+(e«7)<c(e- c ' + e 

£(t) be the smallest integer such that 

1 

c£ £ > lost log log t, 

£ 



(4.23) 



(4.24) 



(4.25) 



(4.26) 



(4.27) 



(4.28) 



(with the condition that 21 + 1 G {2 n — l} n£ N) one sees that (4.27) implies ( 1.15| ). □ 

Appendix A. Some equilibrium estimates 

1.1. A few basic facts on cluster expansion. In this section we rely on the results 
of [9j, but we try to be reasonably self-contained. We let Z 2 be the dual lattice of Z 2 
and we call a bond any segment joining two neighboring sites in Z 2 *. Two sites x, y in 
Z 2 are said to be separated by a bond e if their distance (in M 2 ) from e is 1/2. A pair 
of orthogonal bonds which meet in a site x* £ Z 2 * is said to be a linked pair of bonds 
if both bonds are on the same side of the forty-five degrees line across x* . A contour is 
a sequence eo, . . . , e n of bonds such that: 

(1) ei / ej for every i / j, except possibly when = (0,n) 

(2) for every i, ej and e,+i have a common vertex in Z 2 

(3) if four bonds ei,ej + i and ej,ej+\,i / j,j + 1 intersect at some x* £ Z 2 *, then 
d, ej+i and ej, ej+i are linked pairs of bonds. 
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If eo = e n , the contour is said to be closed, otherwise it is said to be open. Given a 
contour 7, we let A7 be the set of sites in Z 2 such that either their distance (in M 2 ) 
from 7 is 1/2, or their distance from the set of vertices in Z 2 * where two non-linked 
bonds of 7 meet equals 1 / y2. 
We need the following 

Definition A.l. Given V C Z 2 , we let V C M 2 be the union of all closed unit squares 
centered at each site in V , and V be the set of all bonds e G Z 2 * such that at least one 
of the two sites separated by e belongs to V. 

Given a rectangular domain V C Z 2 , a configuration a G fiy and a boundary 
condition r on dV, let a^ T,+ ^ be the spin configuration on Z 2 which coincides with 
a in V, with r on dV and which is + otherwise. One immediately sees that the 
(finite) collection of bonds of Z 2 * which separate neighboring sites x, y G Z 2 such that 

Ox 7^ o-y T ' + ^ splits in a unique way into a finite collection T T (o~) of closed contours. 
It is easy to see that T T (a) n V consists of a certain number of closed contours, plus m 
open contours, where m is such that going along dV one meets 2m changes of sign in r. 
Note that the collection of the 2m endpoints of the open contours is fixed uniquely by 
r. We write ^ T open {o-) for the collection {71, . . . , 7 m } of open contours in T T (a) n V. Of 
course, the open contours ji have to satisfy certain compatibility conditions: 7, and 7^ 
have no bond in common if i j, and if they meet at some x* G Z 2 *, each of the two 
linked pairs of bonds belongs to only one contour. Moreover, each ji is contained in V 
and the collection of the endpoints of the {7i}i< m must coincide with that dictated by 
r. We will write {71, . . . ,7 m } ~ t to indicate that the collection of open contours is 
compatible with r. 

The following result can be easily deduced from [91 Sec. 3.9 and 4.3]. Writing as 
usual TTy for the equilibrium measure in V with b.c. r, one has 

Theorem A.2. There exists /3q such that for every j3 > (3q the following holds. For 
every rectangle V C Z 2 , every b.c. r on dV and every collection {71, . . . ,7 m } of open 
contours compatible with t, one has 




(A.l) 



where the Boltzmann weight ^({71, . . . ,7 m }; V) is defined as 



An(UiA 7i )^0 
7i| is the geometric length of 7$ and 




(A.2) 




(A.3) 



{7l,-,7m}~T 

The potential $ satisfies for every A C V, |A| > 2 and for every x G V : 



$(A)| <exp(-2(/3-/3o)d(A)) 
*({z})|<exp(-8()S-A))) 



(A.4) 
(A.5) 
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where, for connected (in the sense of subgraphs of the graph 1?) A, d(A) is the length of 
the smallest connected set of bonds from A ( cf. Definition A.l) containing all the bonds 
which separate sites in A from sites in A c . If A is not connected then d(A) := +00. 

The fast decay property of $ (with respect to both (3 and d(A)) has the following 
simple consequence: 



Lemma A. 3. |9| Lemma 3.10] There exists (3' depending only on (3q of Theorem A. 2 
such that for [3 > /3' , for every bond e € Z 2 * and for every d > one has 

e -2(/9-A))d(A) < e -2(p-p' )d_ 



E 



(A.6) 



AcZ 2 :eeA 
d(A)>d 



This allows to essentially neglect the interaction between portions of a contour which 
are sufficiently far from each other. 

In order to apply directly results from j9] to obtain the estimates we need, we define 
the canonical ensemble of contours. Let a, b be sites in Z 2 . Then, for any open contour 
7 which has a + (1/2, 1/2), 5+ (1/2,1/2) G Z 2 * as endpoints, in formulas a b (with 
some abuse of language, we will sometimes say that 7 connects a and b), we define the 
probability distribution 



V a ,b(i) ■= {Za,b) 1 exp < 



-2/3| 7 |- Y, *( A )\ = (Za,br 1 *(r,' 



l 2 : 
AnA7^0 



and of course 



(A.7) 



(A.? 



T 1. 
7:a<->0 



Note that we do not require that 7 C V and the sum in VP is now over all (connected) 
sets A C Z 2 . The expectation w.r.t. V a ,b will be denoted by £ a ,b- 

1.1.1. Surface tension and basic properties. Let n be a vector in the unit circle S such 
that n • e 1 > and call <f>^ the angle it forms with e% (of course, —ir/2 < (b^ < 7r/2). 
For iVGN, let 6jv,fj = (A^, J/A^.n) S Z 2 where y^fi = max{?/ G Z : y < A r tan(0^)}. Let 
also := (0,0). Then, it is known [91 Prop. 4.12] that, for (3 large enough, the surface 
tension introduced in (1.2) is given by 

Tp(n) ■= - Jim , v logZ 2)6wf ., (A.9) 



(3d(0, b 



N,n) 



where, if x, y £ M 2 , d(x, y) is their Euclidean distance. To be precise, one has to assume 
that (fifi is bounded away from ±7r/2 uniformly in N, but this will be inessential for us 
since we will always have ^ small. 

One can extract from [9j Sec. 4.8, 4.9 and 4.12] that the surface tension is an analytic 
function of (always assuming that (3 is large enough), and by symmetry one sees 
that it is an even function of <j>^. In (9j Sec. 4.12], sharp estimates on the rate of 
convergence in ( A.9| ) (e.g. (A. 13) below) are given. 
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1.2. Proof of (3.20). The domain El(Ql) which appears in (3.20) is a rectangle with 
height shorter than its base, and the b.c. r is + on the South border and — otherwise. 
Since the event that the unique open contour reaches the height of the South border 



of A is increasing, in order to prove (3.20), by the FKG inequalities we can first of all 



move upwards the North border of El(Ql) until we obtain a square (of side 3L, which 
however here we call just L); we let therefore V := {1, . . . ,L} 2 . Secondly (always by 
FKG) we can change the b.c. r to r' > r by first fixing a 5 > and then establishing 
that t' x = + if x = (xi, x-i) £ dV with X2 < \_5I?-' 2+£ \ , and t' x = — otherwise. 

Given a configuration a £ fiy, let 7 be the unique open contour in T T n (a): of 

course, 7 C V and a± 02, where 01 := (0, |5L 1/2+e J) and a 2 := (L, [5L l / 2+£ \). We 
let h(j) := max{x2 : (xi,X2) € 7} be the maximal height reached by 7, while as usual 
e > is small and fixed. Looking at (A.l) and ( |3.20 ), we see that what we have to 
prove is that for every fixed 5 > one has for every L € N 

^7~'' 1 ~ 1 - ..-'-L- (A 10) 



^1 < e -d* 



E(V,t>) ' H(y,r') 

for some c(/3, 5, e) > 0. We will always assume that j3 is large enough. 

First we upper bound the numerator in (A. 10): with the notations of Section [l.l| (cf. 
in particular (A.7|) and setting for a given contour 7 and a given V C 1? 



one has 



*v(7) := E ^ 

AcZ 2 : An A7^0,AnV c ^0 

< 2 ai>a2 £: aiia2 l {h(7)>25il/2+e} exp($ y (7)) 



(A.ll) 



(A.12) 



< Z aim JV aiA Ml) > 2SL^)j£ aua2 [exp(2cD y (7))] 



where in the first step we simply removed the constraint that 7 C V, which is implicit in 
the requirement 7 ~ t' . It follows directly from [3 Prop. 4.15] that the first square root 
is smaller than exp(— cL 2e ) (note that we are requiring the contour to reach a height 
which exceeds by 5L l l 2+£ the height of its endpoints). On the other hand, from O 
Th. 4.16, in particular Eq. (4.16.6)] and the fast decay properties of $ (in particular 
Lemma A.3| ) it is not difficult to deduce that the second one is upper bounded by 
exp (c(log L) c ) . Moreover, one has |£l Eq. (4.12.3)] that 



Z, 



Ol, Q2 



<cC9)- 



-pr {e!)L 



L 



(A.13) 



where of course Tp(e\) is the surface tension in the horizontal direction and we used 
the fact that d(a±,a2) = L. In conclusion, we have 

M < exp [-f3Tp(ei)L - cL 2e ] . (A. 14) 

Next we observe that, again from (9j Th. 4.16 and Eq. (4.16.7)], 

3(V, r) > exp [-j3rp{e x )L - c(logL) c ] (A.15) 

which together with (A. 14) concludes the proof of (3.20). □ 
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1.3. Proof of Claim 3.10[ In this section, V is the rectangle {(i,j) £ 2 2 : 1 < i < 
L, 1 < j < 4[(2L + 1)^' 2+6 ]} and the b.c. r is defined by r x = - for x G A : = 
{(*,0) 6Z 2 : \i — [L/2\ \ < L 3e } and for x = {x 1 ,x 2 ) G <9V with x 2 > 2 f(2L + l) 1 / 2 ^] ; 
Tx = + otherwise. Moreover, C is the infinite vertical column C = {(xi,X2) G M 2 : 
x\ = LL/2J}. Write Ai + (1,0) (resp. A2) for the left-most (resp. right-most) point 
in A. For every a G fiy there are two open contours in (a): 71 and 72, and we 
establish by convention that 71 is the contour which contains Ai + (1/2, 1/2) as one of 
its endpoints. Two cases can occur (see Figure [6]): 

• either Ai S> A2 and w\ 23- 102, where w\ := (0,2[~(2L + l) 1//2+e ~|) and w 2 ■= 
(L,2\(2L + l)V*+*\) ' 

• or wi 2* Ai and A2 23- w 2 - 





IF 






Y, n_ 







Y, 

K151 



A, 



A 2 



^1 

FF 
A, 



A 2 



Figure 6. The two topologically distinct possibilities: either 71 con- 
nects Ai to A2, or it connects w\ to Ai. The fist case is very unlikely, 



see (A. 18) 



Let C\ (resp. C 2 ) be the vertical column at distance [L e \ to the left (resp. to the 
right) of the column C. Then, one has the 



Lemma A. 4. The probability that appears in Claim \3A0 can be upper bounded as 

n^ +A \r c )<nUr c ), (A.16) 

where 



T := {wi & Ai and 7$ n C{ 



1,2}. 



(A.17) 

Therefore, from Theorem |A.2 we see that to prove Claim |3. 10 it is enough to show 
that 



A/i 



£ {7ll72} ~r *({7i,7 2 }; ^)i {Ai 



71 



S(V,r) 



J(V,r) 



A 2 } < e -cL^ 



and that 



A2 



£{7i,72}~r ^({71, 72>; T/)l {Ai 7i wi} l {7l nC7^0} _ cL . iE 



(A.18) 



(A.19) 



S(V,r) • 

for some positive c = c(f3,e). 

Proof of Lemma A.4\ Since the event T c is increasing, we note first of all that thanks 
to FKG we can enlarge the system from A to V and change the b.c. from (— , +, A) to 
r. Secondly, we observe that the event T implies T. □ 
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1.3.1. Lower bound on H(V,t). We will prove that there exists a positive constant d 
such that for (3 large 

3(V, r) > exp (— /3r / g(ei)(L - c'L 3e )) . (A.20) 

Since we want a lower bound, we are allowed to keep only the configurations {71, 72} ~ t 
such that Wi & Aj and ji does not touch the column Cj, for i = 1,2. Call Qi,i = 1,2 
the set of configurations of 7 j allowed by the above constraints. 
Using the decay properties of one sees that 

H(U,r)>c [ E *( 71 ;F) 

The square is due to the fact that 71 and 72 essentially do not interact because their 
mutual distance is larger than L £ (the residual interaction can be bounded by a constant 
which is absorbed in c). It remains to prove that 

E *(7i; V) > eM-(3MSi)((L/2) - c 'L 3e )) (A.22) 



(A.21) 



for some positive d . This is an immediate consequence of Lemma A. 6 below (applied 
with k = e), together with the fact that d(wi,A{) = L/2 — L 3e + 0(L 2e ), of the fact 
that the angle cj) formed by the segment w\/S.\ and e\ is 0(L _1//2+£ ), and finally of the 
analyticity of the surface tension and its symmetry around e\ . 

1.3.2. Upper bound on A/i . Using rough upper bounds on the number of paths 71 which 
connect Ai and A2 and the decay properties of $ (in particular Lemma A.3), one sees 
that for L large 

cL 3e 



A/i < e 



for some c = c((3, e) > 0, where of course one uses the fact that d(Ai,A2) 
Moreover, Theorem 4.16 of [9] ensures that 



£ *(75 V) < exp(-/3r /3 (e 1 )L + c(logL) c ), 



(A.23) 



2L 3e . 



(A.24) 



■y(ZV: wi*->W2 



which, together with (A.20), concludes the proof of (A. 18). 



1.3.3. Proof of (A. 19). The estimate we wish to prove is very intuitive: if the path 71 
makes a deviation to the right to touch the column C\, it has an excess length, and 
therefore an excess energy, of order L 3e with respect to typical paths. The actual proof 
of (A. 19) is a straightforward (although a bit lengthy) application of results from [9] 
and of the FKG inequalities. We sketch only the main steps. 

First of all, letting ^(71,72) := min{d(xi, X2), Xi G = 1,2}, we show that the 
contribution of the configurations such that c#( 7 i, 72) < L £ is negligible. To this purpose, 
decompose first of all A/2 as A/2 = A/^ + A/^' where 



E 

{71:72}'- 



^({71,72}; U)l {Ai 7i^ } l {7inCl ^ 0} l {(i(7li72)<L£} . (A.25) 
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Consider the paths 7, as oriented from Wi to Aj and, if (£(71, 72) < L £ , call P := 
P (71,72) '■= S Z 2 * x Z 2 * where xi is the first point in 71 n Z 2 * which is at 

distance less than L e from 72, and x 2 is the first point in 72 HZ 2 * at distance less than 
L e from x\. Of course, P can take at most 1? different values (this is a rough upper 
bound) and we can decompose J\f 2 as M 2 = X) P -^2,p where M 2>p contains only the terms 
such that -P(7i,72) = P- Given (71,72) such that ^(71,72) = P, for i = 1,2 one can 
write 7i as the union of 7^ and 7", where 7^ connects tOj to Xj, and 7" connects Xj to 
Aj. Using the decay properties of <& one sees that, uniformly in p and in {7,'}i=i i 2 5 

]T *({ 7 i,72};V) <c*K;y)*(^;n (A.26) 

where the sum runs over all the configurations of {7" }j=i 2 compatible with {7^)1=1 2- 
Let £ be the set of paths 73 which connect x\ to X2, and such that the concatenation 
of 7^,73 and j' 2 1S an admissible open path, call it simply 7, connecting w\ to w% and 
contained in V. Of course, the set £ depends on {t£}i=i 2- Then, one sees that 

2 *({71,72>;F) < e cL£ J] *( 7 ;F). (A.27) 

{T?}i=i,a ^ eS 

In conclusion, summing over the admissible configurations of {7^}i=i2 an d over the 
possible values of p, recalling (A. 24) and the lower bound ( A.20[ ), we have shown that 



^ <e- cL3£ . (A.28) 



E(r,V) 

As for J\f 2 , using the decay properties of the potential $ one sees immediately that, 
since (£(71,72) > L £ , the mutual interaction between the two paths can be bounded by 
a constant, so that 

A^'<c Yl *(7i;^)l {71 nc^0}X E VfaV). (A.29) 

■yiCV: Ai+iwi 12CV: /S.2^>W2 



Recalling (A. 21) one sees therefore that 



K <c—Q—,, (A.30) 



E(t,V) ~ (1-Q) 
where 

_ E {7C y : Ai<2>wi} ^(7;^)l{ 7 nc^0} 

and we are left with the task of proving that Q < exp(— cL 3e ). Note that Q is nothing 
but the equilibrium probability iry^CiCi 7^ 0), where 7 is the unique open contour for 
a system enclosed in V and with boundary conditions f given by f x = + for x = (i, 0) 
with i < [L/2\ - L 3e and x = (0,i) with i < 2[(2L + l) 1/2+£ J, and f x = - otherwise. 
Morally, one would like to apply O Th. 4.15] to say that Q < exp(— cL 3£ ); such 
result however cannot be applied directly because of the entropic repulsion effect that 
7 feels due to the South border of V, and we need to take a small detour. Consider 
the L-shaped domain W obtained as the union of the rectangles V and V, where 
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V' = £ I? : -L 1 ^ 6 < j < 0, 1 < i < [L/2\ - L 3e - 1}, with boundary 

conditions f given by f = f on dW n dV and f' = + on dW n 91/', see Figure [7| 

y 



w +(1/2,1/2) 



(1/2,1/2) 




Figure 7. The L-shaped domain W (for graphical convenience, pro- 
portions are not respected in the drawing) with its boundary conditions 
f'. For the construction of 7', one should imagine that the spins in 
the framed region are set to — . The sites marked by * denote the *- 
connected set A + (7 ; ). The drawn configuration of 7 is entirely above the 
straight line going through w\ + (1/2, 1/2) and Ai + (1/2, 1/2), i.e. the 



spin configuration a belongs to the set T" appearing in (A.36). 



Below we will prove 
Lemma A. 5. One has 



Q = 4( 7 n d + 0) < 4(7 nc^ 0|r') < 



4(7nCi^0) 



(A.32) 



where V = {a € £lw '■ 3 inside V a ^-connected path of + spins which connect the site 
Ai + (0,1) to one of the sites (l,i) imt/t 1 < i < 2[(2L + l) 1 / 2 ^]}, see Figure [?| 



The numerator in the right-hand side of (A.32) is smaller than exp(— cL 3£ ). Indeed, 
it suffices to remark that (cf. the notation ( A.7[ )) it is smaller than 

ZwtAi [ 1 { 7 nci^0}exp($ w (7))] y/V Wl ,a x (7 n Ci / 0)^i,Ai ( ex P ( 2$ w(7))), 



ftoi,Ai [exp ($w(7))] ^i,Ai [exp (#w(7))] 

where $^(7) was defined in ( |A.ll ). Theorem 4.15 of [9] says directly that 

PwiAihnC! /0) <exp(-cL 3e ), 



-(A.33) 
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while the fast decay of $, together with [0, Th. 4.16], implies that 



£ Wl Ai [exp (2$ w (7))] < exp(c(logL) c ) 
SwiAi [ ex P($wK7))] > exp(-c(logL) c ). 



(A.34) 
(A.35) 



Roughly speaking, typical paths (under V W1) a 1 ) have a small intersection with W c 
(again, the precise estimates follow from Th. 4.15]). This is why we enlarged V to 
W: if W were replaced by V, the intersection would not be small any more and the 
expectations in (A.34 )- (|A.35| ) would not be under control. 

The denominator in ( A.32[ ) is also not difficult to deal with: one observes (see Figure 
[7]) that the event V is implied by the event T" ={7 does not go below the straight line 
which goes through A\ + (1/2,1/2) and w\ + (1/2,1/2)} (we will write symbolically 
7 > (Ai^i)). Indeed, the subset of A7 where spins are + is *-connected and satisfies 
the requirements of V . Therefore, 7Ty(T') > exp(— cL £ ). Indeed, 



4'(r')>^(r") 



{7>(Aiiui)} 

£w*(7;W) 



(A.36) 



the numerator is lower bounded by 



expf-^ra^Aj^^i, Ai) - c(d(w 1 , Ai)) £ 



via Lemma A. 6 (take k = e/2) and the denominator is upper bounded by 



exp[-/3T /3 (5 !i , 1 A 1 )d( u, i 5 A i) + c(logd(wi, Ai)) c ] 

via [{JJ Th. 4.16], where v Wi a 1 is the unit vector pointing from w\ to Ai. 

Summarizing, we have obtained Q < exp(— cL 3e ) and, via (A. 30) and (A. 28), we 
have proven ( A.19| ). 



Proof of Lemma A. 5. Given a configuration a £ Qw-, imagine to replace all its spins 
in dV PI V by — , cf. Figure [7J then, associated to the restriction ay G f2y, there are 
exactly two open contours in V. The endpoints of these two contours are (1/2,1/2), 
wi + (1/2,1/2), Ai + (1/2,1,2) and Ai + (-1/2,1/2). Under the assumption that 
a G r', one sees immediately that one of the two contours connects w\ to Ai (this is 



nothing else but the open contour which we have called 7 so far, e.g. in ( A. 32 )); we will 
call 7' the second open contour, see Figure [7| Given a possible configuration for 7', V 
is divided into two components, call them V^(t / ) ) where V~ {^') is the one "in contact 
with" V . It is clear that the intersection A + (7') := A7 7 D V + {~f') is a *-connected set 
{i.e. any two of its points can be linked by a *-connected chain belonging to A + (7')) 
and all spins are + there. It is important to remark that if we take a G V and flip any 
spin in V™ 1 := V + (-y') \ A + (7 ; ), the configuration of 7' does not change. Also, if (with 

abuse of notation) we let iry denote the equilibrium measure in Vy 1 * with b.c. + on 
the portion of the boundary which coincides with A + (7 ; ) and f otherwise, one has 



vry( 7 nCi /0) >4( 7 nCi /0), 



(A.37) 
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by FKG since the event 7 Pi C\ 7^ is increasing. One has then, with S the set of 
possible configurations of 7', 



<(T') 

^7r^(7nC 1 /0|r';7 / =O 



(A.38) 



<(r') 



5> $ ( 7 nCi#0) 



<(r') 



where we used (A.37) in the second inequality. 



□ 



1.3.4. A technical lemma. Let a := (01,02) £ Z 2 * and 6 = (61,62) £ with b\ > a\. 
Let Uafe be the unit vector pointing from a to b and <j> a b be the angle which v a f, forms 
with ex. Assume that — k/4 < <j) ab < it /A. Let A > 0, k > 0, let U a ^ = f/ a ,&(A «) C M? 
be the cigar-shaped region which is delimited by the two curves 



(x - ai)(6i - x) 



ax 



1/2+K 



, x G [ax, bx], 



and t/^Y be the upper half of U a ^, obtained by slicing U a j, along the segment ab. Also, 
we will denote by £ Qi & = T, aj b(A, k) the set of all open contours 7 having a and b 
as endpoints, and such that every bond in 7 has non-empty intersection with U a ,b'i 
similarly we define Then, 

Lemma A. 6. Let (3 be large en ough, and consider a domain V C Z 2 such that V 
contains U~}~,(A,k) (cf. Definition A.l). There exists c depending on (3,A,k such that 



*(7! V) > exp [-PT/3(vab)d(a, b) - c(d(a, b)f K } 



(A.39) 



This result can be obtained via a repeated use of Theorem 4.16 of [9] . The error term 
exp(— c (d(a, b)) 2K ) is very rough (but sufficient for our purposes) and can presumably 
be improved. We do not give full details because they are a bit lengthy, although 
standard, but we sketch the main steps. 

First of all, let for simplicity of notations L := b\ — a\ and A 1 := A/1Q. Then, one 
proceeds as follows (keep in mind Figure |8| : 

• for every — n < i < n, with n = log 2 (L) — 2, let z% = (xi, yi) be a point in Z 2 * 
at minimal distance from {xi,^ b . A , K (xi)), where 



ax + (61 - ax) 



, ■ Ifi 

1 signfo) 2 _ 

4 

j=0 



(A.40) 



• remark via elementary geometrical considerations that for every — n < i < n, 
the cigar-shaped set U Zi:Zi+1 (A' , k) is entirely contained in U^ b (A, k); 
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U x , x (A',k) 

-1 



U„. h (A,k) 



<x,x (A',k) 

-2 -1 




U a , b (A',k) 



Figure 8. A typical path 7 which contributes to the lower bound 
(A. 39). For graphical convenience, we have assumed that a and b 



have the same vertical coordinate, and not all the cigar-shaped sets 
U Zit z i+1 (A' , k) have been drawn. 



restrict the sum (A. 39) to the paths 7 which, when oriented from a to b, go 
through the points z- n , z- n +i, ■ ■ ■ ,z n (in this order), and such that the portion 
of the path between z% and z%+i belongs to Y< ZuZi+1 (A' , «); 
remark that, via the decay properties of the potential <£, the interaction between 
two adjacent portions of 7 just defined can be bounded above by a constant; 
apply Theorem 4.16 of [§] to write that for every — n < i < n one has 

J2 V(.r,V) > exp [-pTp(v ZitZi+1 )d(zi,z i+ i) - c(logd( Zi ,z i+1 )) c ] , (A.41) 



z i+l 



(A',k) 



for some constant c depending on A, k, (3. As for the two portions of 7 from a 
to z- n and from z n to 6, they give a multiplicative contribution of order 1 to 
(A. 39) (this is because d(a, Z- n ) = O(l) and d(b,z n ) = O(l), as is immediately 



seen from the definition of n) ; 
• put together the estimates on the contributions coming from the 2n + 3 por- 
tions of 7 obtained in the previous point: using the convexity and smoothness 
properties of the surface tension r«(-) ; one obtains the claim of the lemma. 
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